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ABSTRACT 

We present the first detection of CS in tlie Antennae galaxies towards tlie 
NGC 4038 nucleus, as well as the first detections of two high-J (5-4 and 7-6) CS 
lines in the center of M 82. The CS(7-6) line in M 82 shows a profile that is 
surprisingly different to those of other low-J CS transitions we observed. This 
implies the presence of a separate, denser and warmer molecular gas component. 
The derived physical properties and the likely location of the CS(7-6) emission 
suggests an association with the supershell in the centre of M 82. 

Subject headings: astrochemistry — ISM: molecules — submillimeter — galaxies: 
individual (Antennae, M 82) — stars: formation 



1. Introduction 



The molecule CS is a good tracer of dense gas (n(H2)> 10^-10^ cm ^) in massive star- 



forming regions in our own Galaxy ( Plume et al.|1992 1997 Andre et al.|2007" ) and in nearby 



galaxies such as the Magellanic Clouds (Nikolic et al. 2007), M 51 and Maffei 2 (Paglione 



et al.||1995D , NGC 253 ( [Mauersberger fc HenkeI||1989||Mauersberger et al.||1989t partin et al. 
2005|, IC 342 and M 82 (iHenkel & Bally||1985l iMauersberger & Henkel||1989t iMauersberger 



et al. 1989). Multi-line studies of CS are crucial for the determination of the average gas 
densities in galaxies since CS transitions have excitation thresholds ranging from 10^-10^ 
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cm-3 for the CS(2-1) line ( |Bronfman et al.||1996D up to ~ 2 x 10^ cm-^ for the CS(7-6) 
transition (Plume et al. 1992). In this letter, we report in the center of M 82 the first detec- 
tions of the CS(5-4) and CS(7-6) transitions and the detection of the CS(5-4) line towards 
the Antennae galaxies (NGC 4038). We have also re-observed lower- J CS lines in M 82. 
These two well-known sources were chosen because interferometric sub millimeter /millimeter 
maps previously obtained have shown high concentrations of molecular gas. High resolution 
^^CO(l-O) maps indicate that the nearby (D = 13.8Mpc, see Saviane et al. 2004) Antennae 
interacting galaxies are likely sites of rapid high mass star formation (Wilson et al. 2000). 
M 82 represents an excellent example of a nearby (D=3.25 Mpc, see Dumke et al. 2001) 



starburst galaxy. Low-J molecular line studies (e.g. Fuente et al. 2006 Seaquist et al. 2006 



Martin et al. 2006) have determined its average gas physical parameters. However, the 



molecular emission from the M 82 nucleus appears to come from multiple gas components. 



In particular the detection of abundant CH3OH (Martin et al. 2006) and HON (Brouillet & 



Schilke 1993 ) indicate the presence of high density gas in the nucleus. Our detection of the 



CS(7-6) line not only clearly confirms the presence of very high density gas but its analysis 
(see Sect. |4]) suggests that it is located in the expanding superbubble likely to be associated 
with supernoave remnants (WeiB et al. 1999 Wills et al. 1999 Yao et al. 2006). 



2. Observations and results 



The observations of the CS(5-4) line in the Antennae galaxies (NGC 4038 nucleus) and 
the CS(7-6) transition in M 82 were performed during the spring 2007 at the James Clerk 
Maxwell Telescope (JCMT). We used a position switching mode under medium weather 
conditions {T225 = 0.16). For observing the CS(7-6) transition = 342.883GHz) towards the 
center of M 82, we used the heterodyne HARP-B multi-beam receiver and the ACSIS digital 
autocorrelation spectrometers with a bandwidth of 250MHz. The receiver noise temperature 
(single sideband mode) for the HARP-B central pixel (corresponding to the center of M 82) 
was ~ 344K. The telescope main beam efficiency and half power beam width (HPBW) under 
this configuration were 0.63 and ~ 14", respectively. The CS(5-4) line (i/= 244.936GHz) has 
been observed in NGC 4038 and in M 82 using the receivers RxA3 and A3 of the JCMT, 
respectively, again with the ACSIS 250MHz bandwidth and DAS backends, respectively. 
The receiver noise temperature (double-side band mode) was 240-450K. The telescope main 
beam efficiency and HPBW were 0.69 and ~ 20", respectively. For both M 82 and NGC 4038, 
the pointing and calibration were performed carefully on planets (Mars and Jupiter) and on 
evolved stars. The pointing error was estimated to be < 2" for NGC 4038 whereas is was 
larger for M 82 (< 8-10"). 
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In the case of M 82, we have also detected the CS(3-2) {u = 146.969GHz), the CS(2- 
l)(zy = 97.980GHz) and the CS(4-3) (z/ = 195.954GHz) hnes with the IRAM-30m in wobbler 
switching mode on June 1997 and on May 2008. As backend we used a 512 x 1 MHz or the 
1024 X 4 MHz filterband. For the CS(3-2) line, the system noise temperature was 300K and 
the pointing error 2-4" whereas the CS(2-1) and CS(4-3) lines show Tgys = 281 K and 1655 
K respectively with moderate pointing errors, due to unstable weather conditions. The main 
beam efficiency and the HPBW of the telescope were 0.69 and ~ 17" at 147 GHz, 0.75 and 
~ 25" at 97 GHz and 0.63 and ~ 13" at 195 GHz. 

All the data have been reduced using either the GILDAS or the Specx packages, remov- 
ing the baseline of individual spectra before averaging them. We fitted Gaussian profiles to 
the resulting spectra (Figs, [ijand [2]). In Table [l] As shown in Fig. [ijand Table 1, the 
most striking result is the narrow line width (only 40 kms~^) of the CS(7-6) line in M 82, as 
compared with those (> 200 kms~^) of the lower-J CS transitions (see Sect. [4|. Note that 
both the CS(4-3) and the CS(5-4) intensities and linewidths are dependent on the baseline 
we removed and thus suffer from large uncertainties, estimated to be up to 20% in the total 
integrated intensity. However their velocity positions seem robust even if significantly shifted 
from the expected systemic Vlsr value of M 82 by a factor of 1.3-1.4. Although this charac- 
teristic could be explained by instabilities in the pointing, large pointing differences would 
produce both a shift in velocity and a narrowing in the linewidth, which are not observed. 
The uncertainties are therefore attributed to the contamination of the observations by the 
M 82 South West lobe emission at a velocity of ~100 kms~^ (see the CS(2-1) spectrum in 
Fig. Q. 



3. The physical properties of the high density gas in M 82 and in NGC 4038 



We derive the following integrated line intensity ratios at the center of M 82 : = 
= 0.88, T52 = ^g|2li) ~ 0.37. We did not compute the ratio tj2 = csf2-r) because 
the CS(7-6) linewidth is very different from those found for the other lines and seems to 
imply a different gas component than the one responsible for the lower-J CS transitions (see 
Sect. |4]). Comparing these ratios with those for NGC 253 (see Martin et al. 2005), we found 
substantial differences. The 1^2 and 1^2 ratios in M 82 are around 3 x the values given for the 
NGC 253-180kms~^ component and about 1.8 x the values given for the NGC 253-280kms~^ 
component. These differences probably arise from different gas excitation conditions in the 



M 82 and NGC 253 nuclei (as already suggested by Bayet et al. 2004, 2006). 



Fig. [3] shows the rotation diagram obtained from our observations of CS in M 82, 
corrected for beam dilution and assuming optically thin emission. As the line proffie of the 
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CS(7-6) transition differs significantly from tliose of tlie otlier lines, and as it is expected 
that the gas in the nucleus of M 82 is not in thermal equilibrium, we favour a fit with two 



components, as indicated by the methanol observations of Martm et al. (2006). However 



due to the large uncertainties of the CS(4-3) and the CS(5-4) lines, it is also possible to make 
an acceptable fit of the data with one component. Although we have made a 2-component 
fit, it is likely that there is a continuous change of the excitation temperature. The rotational 
temperatures for the 2-component fit are 4.3 K and 16.3 K, respectively, while for the single 
component fit it is 12.4 K. The rotational temperatures derived from the rotational diagram 
impose a lower limit to the kinetic temperature of the gas. Thus, the detection of the CS(7- 
6) transition may require a kinetic temperature [Tr) larger than ~12-15K. To estimate the 
opacity of the CS lines, observations of isotopes of CS are needed. Unfortunately, no isotope 
detections have been published so far for this position. The effect of the opacity on the 
rotation diagram is difficult to predict, but it is likely to result in an increase of the derived 
rotational temperatures. 

To constrain the density and the kinetic temperature traced by the detection of CS(7-6) 



line in M 82, we explored LVG models (Goldreich &; Kwan 1974 de Jong et al. 1975). Input 



parameters ranged from 30 to 100 K for Tk, from lO^cm"^ to lO^cm"^ for n(H2) and from 
IxlO^^cm"^ to IxlO^^cm"^ for the CS column density (N(CS)). For each model of the grid, 
we fixed the velocity at the CS(7-6) linewidth value, 40kms~^. In these LVG models, the CS 



collision partner is He and the coUisional rates used are those from Lique & Spielfiedel (2007). 



We obtained a good agreement between the observed and modeled CS(7-6) line intensities for 
a Tk range of 65-70 K, a n(H2) of 1.6xlO'5cm-3 and a N(CS) of 1.6x lO^^cm'^. The source- 
averaged total CS column density (from the observations) and ranges from 1.3 xlO^'^cm"^- 
6.7 xlO^'^cm"^, in agreement with LVG model predictions. Comparing our observations with 



the radiative transfer model of the CS molecule developed by Benedettini et al. (2006), the 



CS(7-6) observed integrated line intensity could only be reproduced for Tk > 50K and n(H2) 
> lO^cm"^. Indeed, the best fit was obtained for Ti^=70K and n(H2)=10^cm~'^ (chemical 



age of 2x10 yr, see Benedettini et al. 2006). This analysis clearly shows that emission from 



high- J (at least 7-6) CS transitions reveal the presence of very high density gas in M 82. 

For the NGC 4038 nucleus, we were not able to perform a similar study as no other 
transitions of CS were observed towards this source. However, we make a rough estimate the 
source-averaged total CS column density from the CS(5-4) line in NGC 4038 to be N(CS) = 



1.8xl0i3cm-2. 



In both sources, the CS column densities are consistent with our model 



predictions of gas undergoing high-mass star formation (Bayet et al. 2008) (see Model 15 
listed in Table 9). For a more detailed analysis, in a forthcoming paper, we shall compare CS 
predicted line profiles from chemical models with the observational line profiles presented in 
this Letter. 
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Discussion 



The angular resolution of the CS(7-6) (beam size of 14") does not allow us to resolve spa- 
tially the emitting region. However we can estimate its location from the radial velocity and 
the velocity width of this line by using the kinematic information provided by interferomet- 



ric maps of other molecules. Brouillet & Schilke (1993) presented a 2" resolution HCN(l-O) 
map of the south west part of M 82. From the channel map of the HCN(l-O) emission cen- 
tered around 213.3 kms^^(±40kms~^), it appears that the CS(7-6) emission may arise from 
the HCN(l-O) clump located at ajaooo = 09^^55^52.2^ and (5j2ooo = 69°40'46.1", within our 



beam. Garcia-Burillo et al. (2001 2002) presented high resolution interferometric HCO and 
H^^CO"'"(1-0) maps of M 82. We re-analyzed the data and found that some of the clumps 
show narrow (< 60 kms~^) emission, similar to that of our CS(7-6) line, confirming the 
likely location of the CS(7-6) emission at ~ (-5", -2") from the center. This is the location 



of the expanding molecular supershell (WeiB et al. 1999 Wills et al. 1999; Yao et al. 2006). 
The high density (> lO^cm"'^) and the high temperature (~ 60-80K) may arise from the 
interaction between the expanding supershell and the ambient gas in M 82. 



In summary, we have presented the first detection of a high density tracer (CS) in the 
Antennae galaxies (NGC 4038 nucleus) as well as the first detections of two high-J transitions 
of CS in the center of M 82. We find that multiple molecular gas components in M 82 are 
necessary to explain the observed line intensities. In particular, while low-J CS lines seem 
to arise from relatively low density gas (~ lO^cm"^), the molecular gas traced by the CS(7- 
6) line must be dense (~ lO^cm"^) and warm (~ 70K) and appears to be associated with 
the expanding supershell in M 82. The high density and temperature may be due to the 
interaction between the expanding supershell and the ambient gas. Similar multi-line studies 



Table 1: Observational parameters. 



Source 


Line 


V 


Tsys 


beam 


/ (T^fe dv) 








rms 






(GHz) 


(K) 


size (") 


(Kkms-i) 


(Kkms-1) 


(Kkms-i) 


(mK) 


(mK) 


M 82 


CS(2-1) 


97.980 


281 


25 


13.3±0.3 


221.9±2.7 


225.7±5.4 


55.5 


5.8 




CS(3-2) 


146.969 


303 


17 


11.2±0.3 


219.7±2.8 


211.3±6.4 


50.1 


4.6 




CS(4-3) 


195.954 


1655 


13 


11.1±1.3 


161.0±9.7 


165.3±23.2 


63.4 


25.8 




CS(5-4) 


244.936 


437 


20 


4.2±0.8 


140.0±38.0 


211.2±44.2 


15.6 


10.8 




CS(7-6) 


342.883 


334 


14 


2.2±0.2 


213.8±1.5 


40.1±2.8 


52.5 


14.1 


NGC 


CS(5-4) 


244.936 


246 


20 


1.7±0.1 


1655.0±4.4 


98.7±8.5 


16.3 


6.7 



4038 
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Fig. 1.— Spectra of the CS(2-1), CS(3-2), CS(4-3), CS(5-4) and CS(7-6) lines from top to 
bottom, respectively, measured towards the center of M 82. The observed position corre- 
sponds to the center of M 82 (ajsooo = 09'^55'"51.9" and 5j2ooo = +69°40'47"). The CS(2-1) 
and the CS(4-3) spectra have been smoothed to a common velocity resolution of 6.1kms^^ 
while the CS(3-2), CS(5-4) and CS(7-6) spectra have a velocity resolution of 8.2kms~^, ~ 
lOkms"^ and 6.8kms~^, respectively. The velocity scale (x axis) is expressed in kms~^ units 
while the temperature scale (y axis) is expressed in main beam temperature units (Tmb), 
converted from the antennae temperature via the main beam efficiencies listed in Sect. [2j 
The solid black lines superimposed onto each spectrum represents the Gaussian fit while the 
vertical black line marks the systemic M 82 Vlsr- 



-7- 





-0.01 



1 

-1 


NC 


1 1 1 1 1 
C 4038 

^ IP 


1 1 1 1 


:;s(5-4) ■ 

y 


■ 1 




if-- 




1 







1500 1600 1700 

Velocity (km/s) 



Fig. 2. — Spectrum of the CS(5-4) line towards NGC 4038. The observed position corre- 
sponds to the NGC 4038 nucleus (ajaooo = 12''01'"52.8" and (5j2ooo = -18°52'05"). The 
spectrum has been smoothed to a velocity resolution of 2.3 kms~^. The solid black line su- 
perimposed on the spectrum represents the Gaussian fit while the vertical black line marks 
the systemic NGC 4038 Vlsr- 
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Fig. 3. — Rotation diagram derived from the CS lines towards M 82. The black lines 
represent the linear regression fits for two gas components while the dashed line corresponds 
to one single gas component. The black solid squares show the data with error bars. These 
error bars usually represent the main errors in rotational diagrams which actually correspond 
to those of the integrated intensities (order of 10-20%). This population diagram has been 
corrected for beam dilution effects assuming a source size of 10" for the center of M 82. This 



value is estimated from the interferometic data presented in Brouillet & Schilke ( 1993 ) and 



Garcia-Burillo et al. (2001). 
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for the Antennae galaxies (NGC 4038/39) are necessary in order to determine the origin of 
its high density molecular gas component traced by the CS(5-4) line. In both sources, the 
CS column densities are consistent with model predictions of gas undergoing high-mass star 



formation (Bayet et al. 2008). 



Similar extended multi-transition multi-molecule studies performed on ultra-luminous 



infra-red galaxies by Greve et al. (2006) and Baan et al. (2008) need to be carried out 



on nearby sources. In fact, Miihle et al. (2007) presented a detailed study of the para- 



H2CO in M 82; however this molecule does not trace the dense gas component. Thus, in 
order to compare the physical and chemical properties of the very dense (extragalactic) gas 
between ULIRGs and nearby sources, CS observations (from 2-1 to 7-6) of ULIRGs should 
be performed. 
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